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Abstract 



We have investigated the ground-state properties of carbon isotopes in the 
framework of the relativistic mean-field (RMF) theory. RMF calculations have 
been performed with the non-linear scalar self-coupling of the a meson using 
an axially symmetric deformed configuration. We have also introduced the 
vector self-coupling of the lo meson for the deformed mean-field calculations. 
The results show that the RMF predictions on radii and deformations are in 
good agreement with the available experimental data. It is shown that several 
carbon isotopes possess a highly deformed shape akin to a superdeformation. 
The single-particle structure of nuclei away from the stability line has been 
discussed with a view to understand the properties near the neutron drip line. 
Predictions of properties of carbon isotopes away from the stability line are 
made. 
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I. INTRODUCTION 



Radioactive beams are being used increasingly to produce nuclei at the limits of the 
stability Nuclei near the neutron and proton drip lines are becoming accessible to 

experiments. Abnormally large reaction cross-sections for unstable nuclei near drip lines 
have been interpreted as existence of a large tail of neutron density in the exterior of nuclei. 
Consequently, the so called halo of particles has been hypothesized for nuclei near drip 
lines. The case of n Li has attracted a widespread attention ||. On one hand, properties 
of nuclei with halos and neutron skins are being synthesized and studied experimentally, 
these properties are providing a test bench to probe various theories and models, on the 
other hand. Description of very light nuclei in terms of a core and a set of valence particles 
hovering around the core are proving to be generally successful J7|. For nuclei with a larger 
number of particles, theories with an average field are being used increasingly. Thus, study 
of nuclei with a large range of isospin puts various theories and interactions to test their 
validity. At the same time, it should become possible to discern the isospin dependence of 
the nuclear force by studying nuclei at the extreme limits. 

Relativistic Mean-Field (RMF) theory is one of the more successful ones in recent times 
for describing nuclei with a large range of isospin P)-l2|. Earlier, the importance of the 
p-meson coupling and hence a proper asymmetry energy was emphasized by Sharma et 



al. 0]. Consequently, the RMF force NL-SH proved to be very successful in describing 
various facets of nuclear structure. Owing to a proper spin-orbit interaction in the RMF 



theory, it became possible to describe the anomalous isotope shifts of Pb nuclei ||14j| . It 



was shown |]15|-[T7|] that ground-state properties of nuclei such as binding energies, charge 
and matter radii and deformation properties of nuclei over a large part of the periodic table 
are very well described within the RMF theory using the non-linear scalar self-coupling of 
cr-meson. The Dirac-Lorentz structure of the nuclear force and the nuclear saturation based 
upon the attractive component due to cr-meson and a repulsive component from tu-meson 
make the RMF theory an attractive tool to study nuclear properties. An inherent spin-orbit 
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interaction arising from the Dirac-Lorentz structure of nucleons gives rise to appropriate 
shell effects. Such shell effects are duly responsible to explain the behaviour of isotope shifts 
in Pb and other nuclear chains fT^JT5|j . Based upon the isospin dependence of the spin-orbit 
interaction in the RMF theory, the Modified Skyrme Ansatz (MSkA) [18|] was proposed. 
Herein, the spin-orbit potential was proposed to contain only one-body contribution. The 
MSkA [pl| was shown to exhibit shell effects which are in accord with the experimental data. 
Consequently, isotope shifts in Pb nuclei could be well reproduced with in the MSkA. 

With the RMF theory having been developed to possess several advantages over the 
non-relativistic theories, we study here the ground-state properties of the chain of carbon 
isotopes. The carbon isotopes away from the stability line have been the focus of experimen- 
tal study whereby interaction cross-sections for heavier carbon isotopes have been measured 
19| . The matter and hence neutron radii of very neutron-rich isotopes have thus been de- 



duced. The nuclear structure of nuclei away from the stability line and especially those in 
the vicinity of drip lines is not yet fully understood. A few attempts have been made to 
explore theoretically the ground-state properties of the carbon isotopes [20]. In the present 
work, we use the RMF theory to investigate systematically the nuclear structure of even-even 
carbon isotopes from proton drip line to the neutron drip line. We employ the non-linear 
scalar potential of the a-meson as well as the model with the non-linear self-coupling of 
both the a and u mesons. We will discuss the nuclear structure of carbon nuclei away from 
the stabiliy line. A comparison will also be made with some results obtained using the 
relativistic Hartree-Bogolieubov approach. It will be interesting to see to what extent the 
mean-field approach will be able to describe the properties of known light nuclei such as 
carbon isotopes. In Section II, we present the formalism of the RMF theory and associated 
models of the a and u meson potentials. The ensuing section provides details of the RMF 
calculations. In section IV we present the results and in the last section we will summarize 
our conclusions. 
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II. RELATIVISTIC MEAN-FIELD THEORY 



The starting point of the RMF theory is a Lagrangian density || where nucleons are 
described as Dirac spinors which interact via the exchange of several mesons. The Lagrangian 
density can be written in the following form: 

£ = $(ip - M)tp + -dpo&o - U{a) 

- g a ^o^ - g^tjjtj} - g P ^t>ri) - e^J^ip. (1) 

The meson fields included are the isoscalar a meson, the isoscalar-vector uj meson and the 
isovector-vector p meson. The latter provides the necessary isospin asymmetry. The bold- 
faced letters indicate the isovector quantities. The model contains also a non-linear scalar 
self-interaction of the a meson : 

U{a) = \ml^ + l -g 2 a* + \g^ (2) 



The scalar potential (2) introduced by Boguta and Bodmer [21] has been found to be useful 
for an appropriate description of surface properties, although several variations of the non- 
linear a and uj fields have recently been proposed J22 ] . We have also included the vector 



self-coupling of the uj- meson introduced by Bodmer ||23|| . The corresponding term in the 
Langrangian is represented by the coupling constant g^. Here M, m CT , m w and m p denote the 
nucleon-, the a-, the to- and the p-meson masses, respectively, while g CT , g^, g p and e 2 /47r = 
1/137 are the corresponding coupling constants for the mesons and the photon. 

The field tensors of the vector mesons and of the electromagnetic field take the following 
form: 

R H» = QVp v _ QVff _ g^p x p) (3) 
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The variational principle gives rise to the equations of motion. Our approach includes the 
time reversal and the charge conservation. The Dirac equation can be written as: 

H«V + F(r) + p + 5(r)]}^= (4) 

where V(r) represents the vector potential: 

1 + r 3 

V(r) = ^o(r) + g p T 3 po(r) + e - A (r) (5) 
and S(t) is the scalar potential: 

S(r) = g<r a(r) (6) 

the latter gives rise to the effective mass as: 

M*(r) = M + S(r) (7) 

The Klein-Gordon equations for the meson fields are time-independent inhomogenous equa- 
tions with the nucleon densities as sources. 

{-A + m 2 a }a(r) = -g a p s {r) - g 2 (r 2 {r) - g 3 o- 3 {r) 
{-A + m 2 Ju (r) = gu,p v (r) ~ ^o( r ) 

(8) 

{-A + m 2 p }p (r) = g p p 3 (r) 
-AA(r) = ep c (r) 

where Uo(r) and po(r) are the time-like components of the uj and p meson fields. The 
corresponding source terms are 



A _ 

8=1 



Pi- = E %■ % 

i=l 



(9) 



Z JV 

P3 = E^Pp % ~ tpn 
p=l n=l 



Z 

Pc= i>p 
p=l 



where the sums are taken over the valence nucleons only. In the the present approach we 
neglect the contributions of negative-energy states (no— sea approximation), i.e. the vacuum 



is not polarized. The Dirac equation is solved using the oscillator expansion method |24| . 

The centre-of-mass correction to the total energy is included by taking the centre-of-mass 

energy based upon the harmonic oscillator prescription as given by 
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E cm = (10) 



III. DETAILS OF CALCULATIONS 

The RMF calculations have been performed in a cylindrical basis where axial symmetry 
has been maintained. The fermionic and bosonic wavefunctions are expanded in a basis of 
harmonic oscillator. For the expansion, we have taken 12 shells for the fermionic as well as 
bosonic wavefunctions. 

The pairing has been included using the BCS formalism. We have used constant pairing 



gaps which are obtained using the prescription of Moller and Nix |25j as given by 

A - 4 - 8 

(11) 

A - 4 - 8 
p ~ 'z T 73 

Only a few empirical data on the pairing gaps are available in this chain. The pairing 
gaps obtained for a few nuclei from the experimental binding energies are found to be 



consistent with the above prescription. The centre-of-mass correction is included by using 



the zero-point energy of a harmonic oscillator as in Ref. ]24] 

We have used the RMF force NL-SH jn| for the Lagrangian with the non-linear scalar 
coupling. This force has been employed widely to describe properties of several chains of 
nuclei. It is known to provide excellent results for nuclei on both the sides of the stability 
line. This force has been found to be especially useful for nuclei far away from the stability 
line. 

For the Lagrangian with the vector self-coupling, we have used the force TM1 |26| . 



We have also used a newly developed force NL-SV1 ^7J which has been obtained by an 
exhaustive study of the ground-state properties of nuclei within the framework of the non- 
linear scalar and non-linear vector self-coupling. The details about NL-SV1 will be provided 
elsewhere. The parameters and coupling constants of the forces we have used are given in 
Table I. 

IV. RESULTS AND DISCUSSION 
A. Binding Energies 

The binding energy of carbon isotopes in the deformed RMF calculations with the non- 
linear scalar force NL-SH and with the non-linear scalar-vector forces TM1 and NL-SV1 
are shown in Table II. The values correspond to the lowest energy for the ground state. 
As carbon isotopes are perceived to be deformed, the RMF minimization was performed 
both for a prolate and an oblate shape. The corresponding deformation for the lowest 
minimum and a secondary minimum (if existent) will be discussed below. The binding 
energy for the secondary minimum are given in parentheses. Experimental values from the 
mass compilation of Audi and Wapstra [2£| are also shown for comparison. 



It is seen that the force NL-SH describes the binding energies of light carbon isotopes 
well. For the heavier isotopes, the force NL-SH gives a slight overbinding. In comparison, 
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The force TM1 gives an equally good description for the lighter nuclei. It is observed that 
NL-SH which is a force with a scalar self-coupling of a- meson only and TM1 which includes 
also the self-coupling of cu-meson produce binding energies which are very close to each other 
for light as well as medium-heavy carbon nuclei. Whilst NL-SH lends an additional binding 
of about 1-2 MeV to the heavier carbon isotopes as compared to the experimental values, 
TM1 provides a stronger overbinding to heavier carbon nuclei as compared to the NL-SH 
results as well as the experimental values. 

The scalar- vector force NL-SV1, on the other hand, underestimates the binding of the 
very light carbon isotopes slightly. For the other carbon isotopes, the agreement of the 
NL-SV1 binding energies with the experimental values is qualitatively better than those of 
NL-SH and TM1. A comparison of predictions of various forces with the experimental data 
shows that the RMF theory provides a good description of the binding energies of carbon 
isotopes. 

B. The Quadrupole Deformation 

The quadrupole deformation fa an d the quadrupole moment Q 2 obtained from the rela- 
tivistic Hartree minimization with various forces are shown in Fig 1. The fa and Q 2 values 
are given in Table III. The deformation parameters for NL-SH (Fig. l.a) show that in the 
lowest-energy state the nucleus 10 C is highly prolate {fa ~ 0.54). This nucleus also exhibits 
a secondary minimum with an oblate shape {fa ~ —0.16). The behaviour of fa using the 
forces with the scalar- vector coupling is very different for this nucleus. Both the forces TM1 
(Fig. l.b) and NL-SV1 (Fig. l.c) predict a well deformed oblate shape with fa ~ —0.29 
and fa ~ —0.21, respectively, for this nucleus. It can be seen from Table III that a highly 
deformed (akin to superdeformation) secondary minimum is also exhibited by this nucleus 
with TM1 {fa ~ 0.64) and NL-SV1 {fa ~ 0.58). 

The stable nucleus 12 C is described as oblate shaped by all the three forces. It is shown to 
be less deformed with NL-SH as compared with TM1 and NL-SV1. TM1, in particular, gives 
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this nucleus a strogly deformed oblate shape (fa ~ —0.39). Whilst TM1 and NL-SV1 give 
a single well-defined oblate minimum for 12 C, NL-SH also predicts a spherical secondary 
minimum only about 100 keV above the lowest minimum (see Table II). A comparison 
of the fa values from RMF with the available experimental information will be in order. 
The quadrupole deformation for 12 C has been estimated using various experimental probes. 



Analysis of inelastic a-scattering experiment |29j led to a value of fa ~ —0.29. However, 



another experiment |3(j with inelastic a-scattering reported a value of fa = —0.40. In 
comparison, a value of fa ~ —0.41 has been deduced [[yj using inelastic scattering of triton 
beam. Similar oblate deformations have also been deduced from inelastic electron scattering 
[ 32 1 . All these experiments demonstrate unequivocally that 12 C has an oblate shape in its 



ground state. The RMF theory with both the scalar self-coupling as well as with scalar- 
vector self-coupling models describes the deformation of 12 C very well. The deformations 
predicted by the scalar-vector coupling models support several experimental deductions. 

The nucleus 14 C with the magic neutron number N=8, and its neighbour 16 C are both 
shown to be spherical with all the three forces. However, for nuclei above N=10, the influence 
of the magic number N=8 diminishes rapidly and consequently the nuclei 18 C, 20 C and 22 C 
take up a well-deformed oblate shape in RMF calculations with all the three forces. It is 
seen that there is an abrupt onset of deformation as 2 neutrons are added to the spherical 
nucleus 16 C. The resulting deformation for 18 C is predicted to be less than that for 20C for 
any given force. The nucleus 20 C with N=14 is well in the middle of the shell where the 
largest value of fa is obtained (see Table III). A further addition of a pair of neutrons also 
brings about a well-deformed oblate shape for 22 C. The magnitude of fa for 22 C is lower 
than that for 20 C. As the magic number N=20 approaches, nuclei assume a spherical shape. 

A comparative look at the values of fa and Q2 amongst the three forces shows that for 
nuclei A=18-22, the force NL-SH produces the lowest deformation, whereas TM1 is shown 
to provide the largest value. The force NL-SV1, on the other hand, gives a value of fa which 
lies between that of NL-SH and TM1. Whilst TM1 and NL-SV1 predict shape coexistence 
only for 10 C, NL-SH is shown to give a secondary minimum with another shape for several 



nuclei. 

A comment on the deformation properties of carbon isotopes in various mean-field the- 
ories will be appropriate to make here. The mean-field theories, in particular, those of 
the density-dependent Skyrme type have generally given a spherical shape to the nucleus 
12 C. This nucleus is widely perceived to be oblate shaped. The RMF theory predicts such 
a shape and especially the scalar-vector model gives a highly deformed oblate quadrupole 
deformation for 12 C. 

C. Neutron-Proton Deformations 

The difference in the quadrupole deformation of the neutron and proton densities are 
shown in Fig. 2. The figure shows that for 10 C, the deformation for neutrons and protons is 
the same with NL-SH, although the (3 2 value is very large (~ 0.54) in the lowest energy state. 
The model with the scalar-vector self-coupling does show a difference in the deformation 
of neutrons and protons. As discussed above, forces TM1 and NL-SV1 with these models 
predict an oblate shape for 12 C. Both these forces also show the difference f3 n — f3 p as positive. 
This implies that the proton field is more deformed than the neutron field. TM1 predicts 
this difference (~ 0.05) to be slightly higher than NL-SV1. 

The f3 n — f3 p is close to zero for 12 C with all the three forces although TM1 and NL-SV1 
predict a stronger oblate shape than NL-SH. The nuclei 14 C and 16 C are spherical and hence 
there is a vanishing f3 n — f3 p . For the nuclei 18 C, 20 C and 22 C, a marked difference in the 
deformations of neutron and proton fluids is seen. All these nuclei are neutron rich and 
predicted to be oblate shaped with significantly large [3 2 values. A large negative value of 
Pn — (3 P for these nuclei means that neutron deformation is considerably larger than the 
corresponding proton deformation. The neutron number for these nuclei is 12, 14 and 16, 
respectively. The neutron deformation (3 n is accentuated by these numbers being in the 
middle of the neutron shell. The onset of the deformation at A=18 corresponds closely to a 
significant occupation of the il w = [202]3/2 + orbital. At the same time, there is a substantial 
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lowering of this level as A=18 (N=12) approaches. Details of single-particle structure will 
be discussed in Section IV. E. 



D. Nuclear Radii and Densities 

The rms radii for neutron, charge and nuclear matter for the carbon nuclei are shown in 
Fig. 3. The values of radii are presented in Tables IV and V. The charge radius is obtained 
by folding the rms proton radius with the finite proton size. Experimental matter radii 
deduced from total reaction cross sections and neutron radii deduced therefrom Jl9] are also 



given in Table IV. The radii shown in the figures correspond to the shape for the lowest 
energy state. 

The charge radius of carbon isotopes shows a roughly constant value for most of the 
nuclei with all the forces. Only for 10 C, there is a significant increase in the charge radius as 
compared to 12 C. This is due to the onset of proton drip line in going to the more neutron- 
deficient isotopes. For isotopes heavier than 12 C, the charge radius shows only a marginal 
increase as pairs of neutrons are added. Such a behaviour is common to most isotopic chains. 

The charge radii predicted by the various RMF forces are given in Table V. Charge radii 
of 12 C and 14 C from elastic electron scattering data |33| are also shown in the table. The 
experimental value of the charge radius is estimated at 2.47±0.02 fm for 12 C and at 2.56±.05 
fm for 14 C. In comparison the force NL-SH predicts the charge radius for 12 C as 2.66 fm 
and both TM1 and NL-SV1 give it at about 2.7 fm. Thus, the RMF theory overestimates 
the charge radius of 12 C slightly. An explicit centre-of-mass correction might improve the 
predicted value. 

The rms neutron radius shows a steady increase from A=10 to A=14 for all the forces 
(see Fig. 3). However, for nuclei above A=14 there is a sudden increase in value. The 
increase in the neutron radius is shown to be phenomenal for 16 C as compared to 14 C in all 
the models. Such an increase in the rms neutron radius by ~0.5 fm is obviated by a major 
shell gap at N=8. The increase in the neutron radius for the isotopes heavier than 16 C is 
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modest with a successive addition of a pair of neutrons. This increase in the neutron rms 
radius increases the neutron skin of the heavier carbon isotopes. It is seen that the neutron 
radius of 24 C is only marginally larger than its lighter neighbour 22 C. The magnitude of the 
increase in the rms neutron radius of 24 C or of any of its neighbours is not so much as to 
characterize this as a neutron halo. It may be noted that 24 C is predicted to be near the 
neutron drip line in all the forces. Thus, it is concluded that the neuron halo is not present 
in the carbon isotopes near the neutron drip line. 

The matter radii as shown in Fig. 3 reflect much the behaviour of the neutron radii. 
It is noteworthy that the features presented by various rms radii for the carbon isotopes 
are similar for both the models with scalar self-coupling as well as with the scalar-vector 
self-coupling. There is also a similarity between the results of TM1 and NL-SV1. 

A comparison of the predictions of RMF matter and neutron radii with the experimental 
values (Table IV) shows that the RMF values describe the data very well. Predictions of all 
the forces are within the error bars. Even for the heavier carbon isotopes such as 16 C and 
18 C the RMF values show a very good agreement with the experimental data. Thus, the 
RMF theory is able to provide reliable predictions of the properties away from the stability 
line. 

In Fig. 3a, we also show neutron radii taken from the spherical Relativistic Hartree- 



Bogoliubov (RHB) calculations |34] using the force NL3. Herein, the finite-range Gogny 
pairing was included. The RHB calculations take into account the effect of states near the 
continuum appropriately. As the force NL3 is known to give results very close to those of 
NL-SH, a comparison of the NL3 results with those of NL-SH is worthwhile. The comparison 
of the RHB neutron radii with the RMF+BCS results shows that the RHB values are very 
similar to the RMF+BCS ones. Especially, for spherical nuclei such as 16 C and 24 C, the 
agreement between our results and those of the RHB is remarkable. This suggests that the 
effect of the continuum on the carbon nuclei near neutron drip line may not be significant. 
It is, however, not clear why RHB underestimates the neutron radius of 12 C. 

The L=0 component of the vector (p v ) and scalar (p s ) neutron density for the carbon 
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isotopes is shown in Fig. 4 for the forces NL-SH, TM1 and NL-SV1. Both the vector and 
scalar densities show a similar behaviour. The scalar densities are slightly smaller than the 
corresponding vector (baryonic) densities. The difference in the vector and scalar densities 
is representative of the relativistic effects in the nuclear structure. 

The increase in the radial extension of the densities is seen clearly as the neutron number 
increases. The increase in the densities of heavier carbon isotopes above 14 C is mostly in the 
exterior of nuclei, thus contributing to an increasingly large neutron skin as more neutrons 
are addded. However, for the highly deformed carbon nuclei a considerable part of densities 
lies in the exterior of the quadrupole (L=2) component. Figure 5 shows the L=2 component 
of the neutron vector density p v obtained with the RMF forces for the three oblate deformed 
nuclei 18 C, 20 C and 22 C. The negative density implies that the neutron matter is missing at 
the poles as compared to a spherical shape. The corresponding matter is accumulated along 
the equatorial plane. 

As shown above the nuclei 18 C, 20 C and 22 C are predicted to be oblate shaped with all 
the three forces. As we have seen, there are slight differences in predictions of the magnitude 
of the quadrupole deformation between various forces. The L=2 component of the density 
is sensitive to these differences in the deformation. The L=2 densities show a negative peak 
which changes its position slightly about 2.5 fm from one nucleus to another. 

The densities given in Fig. 5a for the force NL-SH show that the nucleus 18 C, which has 
p 2 ^ —0.32, has the smallest L=2 component as compared to 20 C which has f3 2 ~ —0.41. 
For 22 C the L=2 peak is smaller than that for 20 C as for 22 C is ~ —0.31 which is smaller 
than that for 20 C. However, the peak for 22 C is bigger than that for 18 C due to a larger 
quadrupole moment for the latter. 

A comparison of the L=2 densities for the various forces shows that except for a slight 
variation in the depths of the peaks due to varying /3 2 values and quadrupole moments, 
the features of the L=2 densities are very similar. This is primarily due to the reason that 
all the three forces considered here predict an overwhelmingly oblate shape for these nuclei 
and that the variation in the deformation in various models is similar from one nucleus to 
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another. 

We show in Fig. 6 the L=0 component of the proton vector density for 8 C, 10 C and 12 C. 
The lighter nuclei 10 C and 8 C are proton rich and it is interesting to see how the proton 
density changes as one approaches the proton drip line. Calculations with the force NL- 
SH (Fig. 6a) show that the proton density in the interior of the nucleus 10 C decreases as 
compared to that in 12 C. The decrease in the central density in interior of 8 C over its heavier 
counterparts is seen to be substantial. This takes place at the expanse of the proton density 
at the surface. However, the spatial extension of the proton density is not so much as to 
characterize it as a proton halo. This is due to the Coulomb barrier which inhibits formation 
of a proton halo. 

The behaviour of the proton density for the forces TM1 and NL-SV1 with the scalar- 
vector self-coupling is not very different from that of NL-SH. Due to slightly different pre- 
dictions of deformation for 12 C and 10 C by these forces as compared to NL-SH (see Table 
III), densities of these two isotopes with scalar- vector forces is only slightly different than 
with NL-SH. The behaviour of the density of 8 C is similar with all the three forces. It may 
be pointed out that the densities of 8 C are presented here only for a qualitative comparison. 

E. Single-Particle Levels 

In order to visualize the single-particle levels contributing to the evolution of the prop- 
erties of carbon isotopes as a function of neutron number, we show the single-particle levels 
for NL-SH and NL-SV1 in Figs. 7 and 8, respectively. The upper panel shows the levels just 
below the continuum. In the lower panel, we show only a few levels which are not so deep ly- 
ing. The levels are identified by the quantum number Q n consistent with the Nilsson scheme 
p5| . The numbers in the parentheses indicate the quantum numbers [Nn z A]. Occupancy in 
the few highest lying Q orbitals is shown in Fig. 9 for NL-SH and NL-SV1. Here we show a 
comparison only between NL-SH and NL-SV1 which represent the scalar self-coupling and 
scalar-vector self-coupling models, respectively. The single-particle properties of the force 
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TM1 are expected to be similar to that of NL-SV1. 

The lower panel of Fig. 7 shows that for the nuclei 10 C and 12 C which are predicted to 
be deformed (see Table III), the level p 3 / 2 splits into [101]3/2~ and [110]l/2~. Th nucleus 
10 C being predicted as highly prolate with NL-SH, the orbital [101]3/2~ lies higher than 
[110]l/2 _ . As one moves to 12 C, the two orbitals cross and the orbital [101]3/2~ is suppressed 
in energy as a consequence of the oblate shape of this nucleus as predicted by NL-SH. These 
two orbitals are degenerate for 14 C and 16 C both of which are spherical. The degeneracy in 
the two orbitals is lifted again for the heavier carbon isotopes 18 C, 20 C and 22 C which are 
all predicted to be oblate shaped. The splitting between the two Q partners amounts to 
about 10 MeV for these nuclei. Interestingly, the [110]l/2~ orbital is lifted up significantly 
as compared to its degenerate (unperturbed) position. 

The behaviour of the orbital [101]l/2~ corresponding to the p3/2 level, is seen to be 
influenced considerably by the presence of deformation. For the deformed nuclei and partic- 
ularly those with an oblate shape, the orbital [101] 1/2" is lowered significantly in energy as 
compared to the position in spherical case. For the prolate shaped nucleus 10 C the opposite 
is true. 

For the force NL-SV1, the behaviour of the three orbitals below —10 MeV is very similar 
to that for NL-SH, as shown in the lower panel of Fig. 8. As both 10 C and 12 C are predicted 
to be oblate, the crossing of [101]3/2~ and [110]l/2~ is not observed. Moreover, the splitting 
between these two levels is proportionate to the f3 2 value —0.21 and —0.33 for 10 C and 12 C, 
respectively, as can be noticed readily from the figure. 

For nuclei up to A=14 (N=8), orbitals corresponding to the levels lsl/2, lp3/2 and 
lpl/2 are filled successively. For nuclei above 14 C, the levels in the next shell start filling. 
The upper panel of Figs. 7 and 8 show the single-particle levels in this shell. These levels 
lie in the neighbourhood of the Fermi surface which is shown by the dotted line. The levels 
include Vt orbitals for ld5/2, 2sl/2 and ld3/2. As the predictions of NL-SH and NL-SV1 
on the deformations for A > 16 are similar, the single-particle structure in both Figs. 7 and 
8 is similar. 
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The three f2 components for the ld5/2 level are shown by the dotted curve. These are 
degenerate for the spherical shape at A=16. It is seen that as the deformed shape evolves 
for A = 18, 20 and 22 on adding neutrons to A = 16, the orbital [202]5/2 + corresponding to 
ld5/2 with the largest f2 value is suppressed. Part of this lowering is also expected from the 
increase in the depth of the potential well for nuclei richer in neutrons. The counterparts 
with lower Q values e.g. [211]3/2 + and [220]l/2 + are, on the other hand, lifted up in energy 
for the deformed isotopes 18 C, 20 C and 22 C. The three f2 components converge for 24 C which 
is close to being spherical. The [200]l/2 + component corresponding to 2sl/2 level is close 
to being degenerate with the ld5/2 level for the nucleus 16 C. The orbital [200]l/2 + follows 
a mild lowering in energy with an increase in A. 

The orbitals [202] 3/2+ and [211] 1/2+ for ld3/2 level also show a significant splitting 
for the deformed nuclei, the larger Q orbital being lowered in energy. The consequence of 
the splitting between various Q orbitals shown in the upper panel of Figs. 7 and 8 is that 
there is a gap of about 3-4 MeV in the energy levels, which persists for A=18-22. 

Figures 7 and 8 show how the Fermi energy is changing with neutron number. For 16 C, 
the Fermi energy decreases rapidly as compared to 14 C. This is due to the onset of the 
next shell for neutrons. Above 16 C, the Fermi energy decreases slowly and approaches the 
vanishingly small value for 24 C. This indicates the onset of the neutron drip line. Such a 
behaviour is similar for NL-SH as well as NL-SV1. 

We show in Fig. 9 the occupation probabilities for some of the orbitals close to the Fermi 
energy. The obitals are arranged in order of increasing energy. Owing to the similarities in 
the predictions of the deformations, the behaviour of the occupation numbers is similar for 
NL-SH and NL-SV1 for the levels which are significant only for the very neutron-rich carbon 
isotopes. The level [202]5/2 + is filled up significantly already at A=18. The neighbouring 
levels [200]l/2 + and [202]3/2 + (see Figs. 7 and 8) also show a considerable occupancy at 
A=18 and are being filled up increasingly up to A=22. The bunch of three levels [202]5/2 + , 
[200]l/2 + and [202]3/2 + is responsible for maintaining a very large oblate deformation for 
carbon isotopes with A = 18, 20 and 22. The levels [211]l/2+ and [211]3/2+ play only a 
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lesser role for these nuclei. The latter is filled up only for A=24 which is predicted to be 
spherical. The population of the orbital [202]3/2 + 24 C is quenched as this levels comes up 
again near the continuum in its unperturbed position for the spherical shape. 

In Fig. 10 we show the shapes of a few highly deformed oblate and prolate configurations 
of carbon isotopes. The vertical axis in each figure is the axis of symmetry. The nucleus 10 C 
has an oblate shape in the ground state as discussed above. This nucleus is also shown to 
give a secondary minimum with a highly prolate shape. Both NL-SV1 and TM1 forces give 
a similar predictions for the deformation. It is instructive to compare the single-particle 
structure of 10 C for the oblate and the prolate shapes. The single-particle structure for 
the two shapes is shown in Fig. 11. As there are only 4 protons in 10 C, two protons go 
inevitably to the lsl/2 level. The other two protons go to the second shell. For the oblate 
shape, the Q n = 3/2 _ orbital is lower in energy and the splitting between the Q 7 * = 3/2~ 
and Q 77 = l/2~ is small. For the prolate shape, the lower fi 71 " = l/2~ orbital is suppressed 
in energy and contains a sizeable fraction of the nucleon occupancy. The splitting between 
the two Q n orbitals is obviously larger for the case of the prolate shape (due to a larger 
value of fa) than the oblate case. The gap betwen the lsl/2 level and the next shell is 
about 17.5 MeV for the prolate shape. The large deformation fa = 0.59 achieved using the 
force NL-SV1 and a larger value fa = 0.64 obtained for TM1 correspond to the so-called 
super deformation of nuclei whereby a 2:1 ratio of the axes are achieved within a harmonic 



oscillator scheme. This also leads to a creation of new shells in the deformation space p5 
The oscillator frequencies u± : u z for the prolate shape are in the ratio 1.76 : 1, which lies in 
the vicinity of the ideal value 2 : 1 for the oscillator potential. Thus, 12 C seems to conform 
to a superdeformed configuration in the secondary minimum. 

The nucleus 20 C which is highly oblate deformed (fa = —0.44) with the force TM1 is 
also shown in Fig. 10. This value of the oblate deformation is close to the 2:3 ratio of the 
axes, whereby the axis of symmetry is accordingly shorter by this ratio. 

In Fig. 11 we compare single-particle levels for 10 C for the lowest energy state with those 
of the second minimum in energy for NL-SV1. As the second minimum state is very close 
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in energy (only ~ 0.5 MeV above the lowest energy), it is instructive to see the difference in 
the single-particle structure of a shape-coexisting highly prolate and an oblate shape. For 
the highly prolate shape (fo ~ 0.59), both the Q components of lp3/2 level, i.e. [110] 1/2 - 
and [101]3/2~ contribute to the prolate deformation, whereby the level [110]l/2 _ plays a 
significant role. For the oblate shape, the orbital [101]3/2~ along with [101]l/2~ of lpl/2 
contribute to the deformation. Here the role of the level [110]l/2~ for the oblate shape 
becomes minimal. It is interesting to see a readjustment of various Q orbitals for two very 
different shapes but with almost the same total energy. 

V. SUMMARY AND CONCLUSIONS 

We have investigated the ground-state properties of carbon isotopes in the framework of 
the relativistic mean-field theory using the self-couplings of a and uj mesons. Calculations 
performed in the axially deformed configuration show that many carbon isotopes except 
those with magic neutron numbers are significantly deformed. The force NL-SH with the 
scalar self-coupling is shown to give a highly deformed prolate shape for 10 C in the lowest 
energy state. However, the forces NL-SV1 and TM1 with the scalar and vector self-couplings 
predict an oblate shape for this nucleus. An oblate shape for 10 C is also shown to occur with 
NL-SH, however, in the secondary minimum. For 12 C, the quadrupole deformation obtained 
with the various RMF forces is in good agreement with the values obtained from various 
experiments. 

For all other isotopes above A=16, predictions on the deformation of nuclei are similar in 
the scalar self-coupling and scalar-vector self-coupling models. Both the models predict an 
oblate shape for the heavier carbon nuclei. The isotopes 18 C, 20 C and 22 C are predicted to 
be well deformed. The relative magnitude of the quadrupole deformations and quadrupole 
moments is largest for TM1 and it is seen that deformations produced by NL-SH are slightly 
smaller than those of NL-SV1. It is observed that NL-SH exhibits the phenomenon of shape 
coexistence for several carbon isotopes whereas such a feature is presented by TM1 and NL- 
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SV1 only for 10 C where this nucleus is predicted to possess a very large prolate deformation 
in the secondary minimum lying only 0.5 MeV above the lowest enery state with NL-SV1 
and 0.9 MeV above the lowest energy state with TM1. 

The relative quadrupole deformations of the neutron and proton mean fields show that 
for lighter carbon isotopes, the neutron and proton fields have almost comparable defor- 
mations in all the models. However, isotopes above A=16 show significant differences in 
the deformations of the neutron and proton mean fields. The quadrupole deformation for 
neutrons are found to be larger than the corresponding proton deformation for the heavier 
carbon isotopes. The magnitude of the difference in the neutron and proton deformations is 
higher with the forces with the scalar-vector self-coupling than with the scalar self-coupling 
alone. 

The rms matter and neutron radii obtained with the scalar and scalar- vector self-coupling 
model in the RMF theory agree well with the experimental values on the matter and neu- 
tron radii deduced from the total reaction cross-sections. The neutron radii for the carbon 
isotopes show a gradual increase with an increase in the neutron number. Such a behaviour 
is predicted by all the forces we have employed. However, the increase in the neutron radius 
and consequently the radial extension of neutrons in space for very neutron-rich carbon nu- 
clei is not so much as to characterize it as a neutron halo. This conclusion is in accord with 
the earlier results obtained on carbon nuclei using spherical relativistic Hartree-Bogoliubov 
calculations ||34j| . A similar statement can be made for carbon nuclei in the vicinity of the 
proton drip line that the proton halo for these nuclei is suppressed. 

The single-particle structure and the occupancy of levels near Fermi surface show that in 
the midst of a large deformation for the nuclei 18 C, 22 C and 22 C, there is a pronounced gap 
in the shell structure. We have also looked into the levels which contribute significantly to 
a large value of deformation in these nuclei. It is shown that the magnitude of deformation 
in some carbon isotopes is akin to a superdeformation. The isotope 20 C in its lowest energy 
state and 10 C in the secondary minimum exemplify this behaviour. Such large deformations 
are produced predominantly in the scalar-vector self-coupling model. These predictions on 
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deformation properties change only a little from the force NL-SV1 to TM1. 
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FIGURES 

FIG. 1. The quadrupole deformation 02 (upper) and the the corresponding quadrupole mo- 
ment Q2 (lower) of the lowest energy ground-state in the RMF theory with (a) NL-SH (b) TM1 
and (c) NL-SV1. The fa and Q2 (mb) values of the shape-coexistent secondary minimum are 
shown by a circle enclosed by a square. The values for 12 C from Ref. || (diamond), Ref. @ 
(triangle up) and from Ref. pl[ (open circle) are also shown for comparison. 



FIG. 2. The difference (3 n — (3 V in the quadrupole deformations of the neutron and proton 
mean-fields. Points with a secondary minimum are shown by a square. 

FIG. 3. The rms neutron, matter and charge radii of C isotopes in the RMF theory obtained 
with the force (a) NL-SH with the scalar self-coupling. The neutron radius from RHB calculations 



of Ref. [34 is also shown, (b) Radii with TM1 and (c) Radii with NL-SV1 with the scalar and 



vector self-coupling. 

FIG. 4. The L=0 component of the neutron vector (baryonic) density p v (upper) and scalar 
density p s (lower) with (a) NL-SH (b) TM1 (c) NL-SV1. 

FIG. 5. The L=2 component of the neutron vector density p v for a few deformed nuclei near 
neutron drip line with (a) NL-SH (b) TM1 and (c) NL-SV1. 

FIG. 6. The L=0 component of the proton vector density p v compared for proton-rich nuclei 
with (a) NL-SH (b) TM1 and (c) NL-SV1. 

FIG. 7. The neutron single-particle levels with the force NL-SH. The lower panel shows the 
deeper lying f2 orbitals, whereas the upper panel shows levels in the vicinity of the Fermi surface. 

FIG. 8. The neutron single-particle levels with the force NL-SV1. The details are the same as 
for Fig. 7. 

FIG. 9. The occupation probabilities for the highest five levels in the vicinity of the Fermi 
surface, for (a) NL-SH and (b) NL-SV1. 
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FIG. 10. The shapes of a few strongly deformed carbon isotopes. The upper panel shows the 
shape for the lowest energy state of 10 C and 20 C. In the lower panel the prolate deformation for 
the secondary minimum for 10 C is shown. The vertical axis represents the axis of symmetry. 

FIG. 11. The single-particle structure of the lowest energy state (oblate) and the highly 
deformed secondary minimum (prolate) for 10 C using the force NL-SV1. The particle occupation 
numbers are given in the parentheses. 
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TABLES 

TABLE I. The Lagrangian parameters of the forces NL-SH, TM1 and NL-SV1 used in the 
RMF calculations. 





NL-SH 


TM1 


NL-SV1 


M 


939.0 


938.0 


939.0 


77? 


526 05921 


511.198 


510 03488 


ryi , 


783.0 


783.0 


783. 


m p 


763.0 


770.0 


763. 


9a 


10.44355 


10.0289 


10.12479 




12.9451 


12.6139 


12.72661 


9p 


4.38281 


4.6322 


4.49197 


92 


-6.90992 


-7.2325 


-9.24058 


93 


-15.83373 


0.6183 


-15.388 


94 


0.0 


71.5075 


41.01023 
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TABLE II. The binding energies (in MeV) of even-even C isotopes obtained for the lowest 
energy state with the forces NL-SH, TM1 and NL-SV1. The empirical values (expt.) available are 
also shown for comparison. The numbers in the parantheses indicate the existence of a secondary 
minimum in the vicinity of the lowest-energy ground-state. 



A 


ATT QTT 


Ti\/ri 


t\tt cnri 
IN Ij-o V 1 


expt. 


10 


-60.4 (-59.8) 


-60.2 (-59.3) 


-57.9 (-57.4) 


-60.3 


12 


-89.6 (-89.5) 


-90.1 


-87.6 


-92.2 


14 


-106.6 


-106.9 


-104.4 


-105.3 


16 


-112.3 


-112.4 


-109.6 


-110.8 


18 


-117.0 (-116.9) 


-118.6 


-114.4 


-115.7 


20 


-121.8 


-123.4 


-118.5 


-119.2 


22 


-122.7 (-122.4) 


-123.9 


-118.6 


-120.3 


24 


-122.3 


-122.7 


-116.6 
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TABLE III. The quadrupole deformation and the quadrupole moment Q2 (mb) obtained 
in the RMF theory for the C isotopes. 









(h 




Q2 




A 


NL-SH 


TM1 


NL-SV1 


NL-SH 


TM1 


NL-SV1 


10 


0.536 (-0.161) 


-0.294 (0.643) 


-0.213 (0.584) 


85.5 (-25.7) 


-46.9 (102.6) 


-33.9 (93.1) 


12 


-0.23 (0.005) 


-0.388 


-0.328 


-50.6 (0.9) 


-83.9 


-70.9 


14 


0.000 


0.011 


0.005 


0.14 


0.3 


0.14 


16 


-0.005 


-0.006 


0.004 


-1.6 


-2.2 


1.2 


18 


-0.316 (0.387) 


-0.354 


-0.325 


-134.3 (164.3) 


-150.4 


-137.9 


20 


-0.405 


-0.444 


-0.418 


-205.3 


-224.7 


-211.6 


22 


-0.308 (0.029) 


-0.367 


-0.336 


-183.2 (17.1) 


-217.9 


-199.3 


24 


-0.006 


-0.16 


-0.07 


-4.0 


-112.8 


-49.4 
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TABLE IV. The rms neutron radii r n and matter radii r m (in fm) as obtained in the RMF the- 



ory for various RMF forces. The experimental values [19] deduced from total reaction cross-sections 
are also shown for comparison 







r n 














A 


NL-SH 


TM1 


NL-SV1 


expt. 


NL-SH 


TM1 


NL-SV1 


expt. 


10 


2.49 


2.49 


2.48 




2.69 


2.69 


2.67 




12 


2.50 


2.57 


2.56 


2.49 ±0.16 


2.52 


2.59 


2.57 


2.48 ± 0.08 


14 


2.57 


2.59 


2.60 


2.70 ±0.10 


2.54 


2.56 


2.56 


2.62 ±0.06 


16 


3.01 


3.03 


3.04 


2.89 ±0.09 


2.83 


2.85 


2.85 


2.76 ± 0.06 


18 


3.18 


3.22 


3.24 


3.06 ± 0.29 


2.98 


3.01 


3.02 


2.90 ±0.19 


20 


3.32 


3.37 


3.35 




3.11 


3.16 


3.13 




22 


3.47 


3.52 


3.54 




3.25 


3.30 


3.31 




24 


3.59 


3.63 


3.66 




3.35 


3.40 


3.42 
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TABLE V. The charge radius r c obtained with various RMF forces. The experimental values 
for 12 C and 14 C from the electron scattering data |?3| are also shown. 





A 


NL-SH 


TM1 


NL-SV1 


expt. 


10 


2.93 


2.92 


2.91 






2.00 


1.(6 


O 71 
I. 11 


o a i _i_ n no 
lAi ± U.Uz 


14 


2.62 


2.64 


2.63 


2.56 ±0.05 


16 


2.62 


2.64 


2.63 




18 


2.64 


2.67 


2.66 




20 


2.66 


2.7 


2.68 




22 


2.67 


2.72 


2.7 




24 


2.66 


2.72 


2.69 
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